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We propose a simple and versatile model to understand the deviations from the well-known 
Kolmogorov-Johnson-Mehl-Avrami kinetics theory found in metal recrystallization and amorphous 
semiconductor crystallization. We analyze the kinetics of the transformation and the grain size dis- 
tribution of the product material, finding a good overall agreement between our model and available 
experimental data. The information so obtained could help to relate the mentioned experimental 
deviations due to preexisting anisotropy along some regions, to certain degree of crystallinity of the 
amorphous phases during deposition, or more generally to impurities or roughness of the substrate. 



The interest in thin film transistors made of polycrys- 
talline silicon and silicon-germanium has been driven by 
the technological development of active matrix-addressed 
flat-panel displays ||l| and thin film solar cells |^ . In this 
context, the capability to engineer the size and geometry 
of grains becomes crucial to design materials with the 
required properties. Crystallization of these materials 
takes place by nucleation and growth mechanisms: Nu- 
cleation starts with the appearance of small atom clusters 
{embryos). At a certain fixed temperature, embryos with 
sizes greater than a critical one become stable nuclei; oth- 
erwise, they shrink and eventually they vanish. Such a 
critical radius arises from the competition between sur- 
face tension and free energy density difference between 
amorphous and crystalline phases (which favours the in- 
creasing of grain volume) yielding an energy barrier that 
has to be overcome to build up a critical nucleus. Surviv- 
ing nuclei grow by incorporation of neighboring atoms, 
yielding a moving boundary with temperature depen- 
dent velocity that gradually covers the untransformed 
phase. Growth ceases when growing grains impinge upon 
each other, forming a grain boundary. The final product 
consists of regions separated by grain boundaries. This 
simple picture has, however, two problems: On the one 
hand, this theory of nucleation and growth predicts an 
energy barrier far from the experimental value so nucle- 
ation would hardly be probable at available annealing 
temperatures ^ . On the other hand, it is known that in 
crystallization of Si over Si02 substrates, nucleation de- 
velops in the Si/Si02 interface due to inhomogeneities or 
impurities that catalyze the transformation Q. There- 
fore, a theory of homogeneous nucleation and growth is 
not entirely applicable to the referred experiments. 

The transformation kinetics is also problematic. It is 
generally accepted that the fraction of transformed mate- 
rial during crystallization, X{t), obeys the Kolmogorov- 
Johnson-Mehl-Avrami (KJMA) model according to 
which X{t) = 1 — exp(— At™), where A is a nucleation- 
and growth-rate dependent constant and m is an expo- 



nent characteristic of the experimental conditions. Two 
well-defined limits have been extensively discussed in the 
literature: When all the nuclei are present and begin to 
grow at the beginning of the transformation, the KJMA 
exponent, m, is equal to 2 (in two dimensions), and the 
nucleation is termed site saturation. The product mi- 
crostructure is tesselated by the so-called Voronoi poly- 
gons (or Wigner-Seitz cells). On the contrary, when new 
nuclei appear at every step of the transformation, m = 3 
and the process is named continuous nucleation. Plots of 
log[— log(l — X)] against log(t) should be straight lines of 
slope m, called KJMA plots. The vahdity of the KJMA 
theory has been questioned in the last few years 1^, and 
subsequently several papers have been devoted to check 
it in different ways However, those theoretical re- 

sults still leave some open questions: For example, an 
exponent between 2 and 3 is experimentally obtained in 
two dimensions, the KJMA plots from experimental data 
do not fit a straight line in some cases |0 , and the con- 
nection between geometrical properties (grain size distri- 
butions) and the KJMA exponent is not clear. 

In this letter, we show that these questions may be 
answered by assuming that nucleation is heterogeneous, 
not in a phenomenological way as in other proposed mod- 
els Q , but sticking to the basic ideas due to Cahn 14 
and Beck : The material is not perfectly homogeneous 
but contains regions with some extra energy (regions with 
some order produced during deposition, or substrate im- 
purities) at which nucleation is more probable. Accord- 
ingly, we introduce a computational model consisting of 
several simple irreversible rules, with the additional ad- 
vantage that it describes simultaneously space and time 
evolution. Furthermore, it allows us to average over a 
large number of realizations in very short computational 
times as compared to other computer models (see the re- 
cent review by Rollett [ p^ for an overview of simulation 
models of recrystallization). 

The model is defined on a two-dimensional lattice 
(square and triangular lattices were employed) with pe- 
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riodic boundary conditions. Every lattice site (or node) 
X belongs to a certain grain or state, g(x, t) = 0, 1, 2,. . . , 
the state being that of an untransformed region. The 
lattice spacing is a typical length scale related to the 
available experimental resolution. Following the idea 
that the amorphous phase has random regions in which 
nucleation is favored, we choose a fraction c of the total 
lattice sites to be able to nucleate. We term these en- 
ergetically favorable sites potential nuclei. These poten- 
tial sites may be interpreted as random sites on a region 
where order is present, not just an isolated critical clus- 
ter. Initially ^(x, 0) = for all lattice sites x and the 
system evolves by parallel updating according to the fol- 
lowing rules: i) A transformed site remains in the same 
state [q{x,t + At) = q{x,t) ^ 0]; ii) An untransformed 
potential site may become a new non-existing state (i.e., 
crystallizes) with probability n (nucleation probability), 
if and only if there are no transformed nearest neighbors 
around it; iii) An untransformed site (including potential 
sites) transforms to a already existing transformed state 
with probability g (growth probability), if and only if 
there is at least one transformed site on its neighborhood. 
The new state is randomly chosen among the neighboring 
grain states. 




FIG. 1. Computer simulation obtained for two stages of the 
transformation process on a 250 x 250 triangular lattice with 
(a) c = 1 (homogeneous nucleation), n = 0.001 and g — 0.8 
(total run time, 1 s.); (b) c = 0.1, n = 0.001 and g = 0.8 
(total run time, 1 s.). 

For the model parameters, we expect a functional form 
n ~ Q-E„/kBT g^j^j g ^ ^-Eg/ksT ^ where i?„ and Eg 

are the energy barriers of nucleation and growth respec- 
tively. Hence, temperature is implicit in the definition of 
n and g. Figure ^ shows the microstructure at two dif- 
ferent stages for two different sets of parameters. As we 
are interested in this paper in how different nucleation 
conditions yield different KJMA exponents and different 
microstructures for isothermal experiments, we define a 
characteristic time t as the time that a grain needs to 



increase its size by one lattice site, and consequently we 
can put 5=1. The simulation time step is therefore this 
characteristic time r. 

We have simulated 1000 x 1000 triangular and square 
lattices and averaged the outcome of 50 different realiza- 
tions for each choice of parameters (characteristic sim- 
ulation times are about 15 to 45 minutes in a Pentium 
II personal computer). The main results are the follow- 
ing: If c ^ 1, then most sites are potential sites, so new 
grains are able to nucleate at every stage of the transfor- 
mation (continuous nucleation). On the contrary, when 
c ^ 1, and n < \, every potential site nucleates at the 
early stages of the process (site saturation). Obviously, 
intermediate values yield a mixed behaviour. Interest- 
ingly, the model parameters tune the KJMA exponent 
between 2 and 3. It is important to note that for small 
values of c, which would in principle mean that growth is 
by site saturation, low values of n (large energy barriers 
for nucleation) lead to m ~ 3, as in continuous nucle- 
ation. 
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FIG. 2. KJMA plots for different sets of parameters. The 
dotted lines represent the theoretical slopes 2 and 3. 

Other forms of experimental behavior lead to the oc- 
currence of non-straight KJMA plots. We argue that 
this fact may be due, on the one hand, to the decay of 
the nucleation rate when n ^ 1, because some poten- 
tial sites are overlapped by already growing grains; and 
on the other hand, when the potential site concentra- 
tion is c 1, the grains grow independently for times 
lower than a characteristic impingement time, propor- 
tional to the mean grain distance 1/c^/^. Figure || shows 
this fact for several choices of parameters n and c. Note 
that when n ^ 1, the potential sites nucleate during the 
earlier stages of the transformation, so the mentioned 
overlapping of potential sites cannot be the cause of the 
bending of the KJMA plots. Therefore, we must conclude 
that heterogeneous nucleation is not the unique cause of 
the unexpected bending of the KJMA plots, as m may be 
affected by anisotropics or preferential crystalline direc- 
tions yielding growth or nucleation rates that may change 
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locally throughout the material. This agrees with the fact 
m is not a reliable guide to characterize the morphology 
of the evolving grains . 

As we have pointed out, our model provides informa- 
tion about microstructure, i.e., number of grains, mean 
grain area, grain size distribution, and so on. For site 
saturation, Weire et al. proposed a phenomenological 
expression for grain size distributions P{A') ^ 

{A')°'-^a"e-°"^' /T{a), where a « 3^65 and A' = A/A 
is the reduced area. The mean area A changes from one 
process to another, but the normalized distribution is 
the same for all. Analogously, in the case of continuous 
nucleation, a simple expression has been proposed [ pO| : 
P{A') = . Figure ^ shows the good agreement be- 
tween the simulations of our model and these theoreti- 
cal predictions. For intermediate-ranging parameters, a 
continuous evolution is obtained from site saturation to 
continuous nucleation grain size distributions. We thus 
have two elements of comparison between our model and 
experimental results: the KJMA exponent, to, and the 
grain size distribution P{A'). 
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FIG. 3. Histograms of the grain size distribution for differ- 
ent sets of parameters. 

In conclusion, we have presented a simple lattice model 
for crystallization which sheds light on the possible causes 
of the experimental deviations from the KJMA theory. 
Thus, preexisting inhomogeneities in the initial state, 
such as regions with a lesser degree of disorder or im- 
purities, dramatically change the product structure and 
the time development of the crystalline phase. One of the 
remarkable points of our model is its versatility, so other 
ingredients can be simply added to the model rules. We 
postpone the detailed study of heterogeneous growth or 
preferential directions to further research. The main con- 
clusion of this work is that the KJMA exponent is not 
enough to understand and to characterize the crystal- 
lization mode in a specific experiment: Indeed, we have 
shown that conditions close to site saturation and con- 
tinuous nucleation give rise to very similar values of to. 



Therefore, studies of the grain size distribution are indis- 
pensable to identify correctly the crystallization mode. 
We stress that the model rules are physically meaningful 
(alternative proposals can be found in Ref. |2^, but are 
far from being physical because they depend strongly on 
the lattice geometry and the site interactions), and lead 
to experimentally verifiable predictions. Due to its ver- 
satility and short simulation times, it is an easy to repro- 
duce, good, and nonexpensive testbed for the design of 
materials and structures with tailored grain size or shape 
properties. 
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